Experimental demonstration of the stability of Berry's phase for a spin- 1/2 particle 



S. Filipp,i 2.3 Bj- Klepp, 



'■^ Y. Hasegawa,' C. Plonka-Spehr,^ U. Schmidt,'* P. Geltenbort,^ and H. Rauch' 

Austria 



00 

o 

O 
(N 

o 

Q 

0^ 



Oh 



> 

m 

oo 
O 



X 



Atominstitut der Osterreichischen Universitdten, Stadionallee 2, A-1020 Vienna, 
^Department of Physics, ETH Zurich, CH-8093 Ziirich, Switzerland 
^Institut Laue Langevin, Botte Postale 156, F-38042 Grenoble Cedex 9, France 
^ Physikalisches Institut, Philosophenweg 12, 69120 Heidelberg, Germany 

The geometric phase has been proposed as a candidate for noise resilient coherent manipulation of fragile 
quantum systems. Since it is determined only by the path of the quantum state, the presence of noise fluctuations 
affects the geometric phase in a different way than the dynamical phase. We have experimentally tested the 
robustness of Berry's geometric phase for spin-1/2 particles in a cyclically varying magnetic field. Using trapped 
polarized ultra-cold neutrons it is demonstrated that the geometric phase contributions to dephasing due to 
adiabatic field fluctuations vanish for long evolution times. 
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The rapidly increasing capability to control and measure 
quantum states on a single particle level (see e. g. fH] and ref- 
erences therein) demands for decoherence free systems and 
robust manipulation techniques. In order to enable high- 
precision quantum measurements ['5] or quantum information 
processing [3] with long coherence times, a system should be 
de-coupled from the environment except for precisely control- 
lable interactions. As part of the quest for reliable quantum 
gates the geometric phase has attracted renewed attention due 
to its potential resilience against noise perturbations. 

In short, the adiabatic evolution of a quantum system re- 
turning after some time to its initial state gives rise to an addi- 
tional phase factor, termed Berry's phase The peculiarity 
of this phase lies in the fact that its magnitude is not deter- 
mined by the dynamics of the system, i. e. neither by en- 
ergy nor by evolution time, but purely by the evolution path 
from the initial to the final state. A vast number of experi- 
ments have verified its characteristics in various systems isj. 
Several extensions, for instance to non-adiabatic, non-cyclic, 
non-unitaryor non-abelian geometric phases have been inves- 
tigated iSHSIitl- For closed quantum systems the geomet- 
ric phase is theoretically well understood and experimentally 
verified. However, for open quantum systems the situation 
is different in that no general framework has found approval 
yet. Concepts of geometric phases for mixed state evolutions 
have been introduced theoretically | T^, U , 12 1 and inspected 
experimentally 1 13 , 14 , 15. 3 . Also dephasing induced by 



the geometric phase has been studied theoretically for several 
settings Potential advantages of geometric 



quantum gates for quantum information processing have been 
topic of recent investigation 12 li l22l 12311 . Furthermore, high 
fidelity geometric gates are cuiTently used in ion traps [24] 
suggesting Berry's adiabatic geometric phase as favourable 
choice for quantum state manipulations. In |25] it is shown 
that the contribution of the geometric phase to dephasing are 
path-dependent like the geometric phase itself, as experimen- 
tally demonstrated in 12611 . and that they diminish for long evo- 
lution times. 

In this letter we consider the situation of an adiabatic evo- 



lution of a spin-1/2 system and explicitly test the influence of 
slow fluctuations onto the Berry phase. We complement the 
result in [26] by analyzing the influence of evolution time on 
the geometric dephasing using a dedicated ultra-cold neutrons 
apparatus We show that the Berry phase is robust against 
adiabatic fluctuations in the driving field, when the evolution 
time is longer than the typical noise correlation time. 

Consider a spin-1/2 particle exposed to slowly varying 
magnetic fields. The Hamiltonian 



H{t) = -ija-B{t) = -ij(a-Bo{t)+a,K{t) 



(1) 



describes the coupling of a particle to the magnetic field B{t) 
by its spin magnetic moment /j. The magnetic field has mag- 
nitude B{t)= \B{t) I and its direction points along the unit vec- 
tor n{t) = (cos'd(f), sin-d(f)sin(p(f), sin-&(f)cos(p(f))'^. o ~ 
(o,T,Oy,0^)^ denote the Pauli matrices and K{t) stands for an 
additional fluctuating magnetic field along the x— axis. Let 
|i±(f)) denote the time-dependent spin eigenstates of H{t). If 
the system is initially in an eigenstate |i±(0)), it will stay in 
an eigenstate during an adiabatic evolution of the B-field. In 
other words, the B-field direction and the polarization vec- 
tor of the particles' spin s{t) coincide for all times, where 
s{t) = Tr[o • p(f)] for the system being in the state described 
by the density matrix p(f). In spherical coordinates s{t) can 
be written as 

s{t) = i(cos9(f),sin9(f)sin(|)(f),sin0(f)cos(|)(f))'''. (2) 

Its length i = |?| represents the degree of polarization. For a 
pure state p(f) — \s±{t)) {s±{t)\ we have s = 1, but in general, 
interactions with the environment lead to mixed states with 
reduced i < 1 as discussed below. 

Within the adiabatic approximation, where 9(f) « 'd(f) and 
(|)(f) w (p(0^ ^ cyclic variation of the B-field coordinates 'd(f) 
and (p(r) leads only to a change of the phase of the initial 
state |s-i-(0)) after the evolution. The final state \s±{T)) — 
g±i(<t's+<t'rf)|jj_(o)) comprises a dynamical ((^d) and a geomet- 
ric phase, i^d — ~ fo E {t)dt /his determined by the inte- 
grated instantaneous energies E{t) = —/j {s±{t)\a ■ B{t)\s±{t)) . 
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It depends explicitly on the dynamics of the state transport. In 
contrast, the geometric phase (^g is determined by a surface 
integral proportional to the solid angle Q. enclosed by the path 
of the state: (^g = —D./2 for a spin- 1/2 particle. It is there- 
fore independent of energy and time. In particular, the B-field 
we have used in our experiment traces out a path with con- 
stant polar angle and varying (p(f) e [0,27t] .Without fluctu- 
ations {K{t) = 0) the geometric phase evaluates in this case to 
(|)^ = — 7t(l — cos-d). 

Field fluctuations in x-direction during the evolution are 
represented by the term OxK{t) (with equivalent results hold- 
ing also for isotropic noise involving a^ y terms 12511 ). Fluc- 
tuations in the Larmor frequency co^ = l/jB/h are then given 
by 2/jA'(f ) /h, which denotes a Gaussian and Markovian noise 
process with intensity Op and correlation time 1 /F, i. e. noise 
bandwidth F. We assume an upper cut-off frequency of the 
noise F^ax ^ <J3l such that the field fluctuations are adiabatic 
with respect to the Larmor frequency. Later, in the experi- 
ment, this is achieved by adding adequately designed noise to 
the field. Consequently, the variations in the path, and there- 
fore in £2, lead to variations in the geometric phase. The ran- 
dom geometric phase (^g is Gaussian distributed with mean 
value equal to the noise-free case, {(^g) = (|)^. Its variance 



a^(r) depends on the evolution time T and is given by 1125 1 
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The dependence on the factor sin has been tested in fit 
A further intriguing property is that for evolutions, which are 
slow relative to the noise fluctuations (T ^ 1 /F), the variance 
of the geometric phase drops to zero for long evolution times 
as the expression in Eq. ^ reduces to a^^{T) °^ j 13311 . This 
contrasts the behaviour of the variance of the dynamical phase 
that increases linearly in time. 

In our experiment we have used neutrons as a precisely 
manipulable spin- 1/2 quantum system. Exposure to a mag- 
netic field leads to Zeeman energy splitting of 2/j„B with 
/j„ = —9.66 X 1Q~^^ JT^'. The experimental setup is shown 
in Fig. [T] and more details can be found in i27il . Neutrons are 
guided from the ultra-cold neutron source at the ILL high flux 
reactor through magnetized Fe polarization foils, which give 
a degree of polarization of about 90%, to the storage bottle 
(filling stage). Their low kinetic energy prohibits penetration 
through the walls of the bottle. During storage the spin orien- 
tation of the dilute (w 1 neutron/cm^) gas of practically non- 
interacting spin- 1/2 particles can be arbitrarily manipulated 
by magnetic fields produced by a 3D Helmholtz-coil setup 
(manipulation stage). The resulting spin polarization is sub- 
sequently analyzed by a combination of a fast adiabatic 7t— 
flipper (ss 99% efficiency) and the polarization foils before 
hitting the detector (emptying stage). A full storage cycle of 
filling, manipulating and emptying lasts about 70 seconds. 

In the following we focus on the 'manipulation stage': 
First, we have measured the adiabaticity of the transport of 
the polarization vector. We have compared the initial to the fi- 
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FIG. 1: Experimental setup: Neutrons passing through polarization 
foils are guided to the storage bottle, where they are stored for typi- 
cally 10 seconds before letting them fall down to the detector. During 
storage their spin is manipulated by magnetic fields produced by the 
Helmholtz coils. 



nal polarization after a cyclic variation of the magnetic guide 
field (B(0) = B{T)) with the neutron spin initially aligned with 
the static magnetic guide field in the negative z— direction. In 
order to find identical initial and final polarization the adia- 
baticity condition requires the rate of change of the B-field 
\{dB/dt)/B\ to be much smaller than the Larmor frequency 
COl of the system. Within the accuracy of the experiment this 
applies for a typical rate of change less than approximately 
0.2(Oi setting an upper limit for the following measurements. 

Secondly, a Ramsey-type interferometric scheme similar to 
has been employed for the measurement of Berry's phase. 
The pulse sequence of the B-field during the manipution stage 
is shown in Fig.|2] The actual evolution is preceded by a 7t/2- 
rotation induced by a rf-field in x— direction with amplitude 
1.6 /jT, duration 10.7 ms and a frequency resonant with the 
magnetic guide field (0) = — 1 /jT. Starting from the eigen- 
state |5+(0)) this generates an equal superposition of spin- 
up and spin-down states, |v(/(0)) = (|i+(0)) + |i_(0))) /%/2. 
A subsequent adiabatic and cyclic B-field evolution of dura- 
tion T induces a relative phase between the states of <^{T) — 
(^g + The resulting spin polarization s{T) can be ana- 
lyzed by 7t/2-pulses, which are offset in phase by zero or 
90° relative to the preparatory 7t/2-pulse, thus yielding a ro- 
tation of s{T) about the x— or y— axis, respectively. A fur- 
ther 7t-flip can be induced with high efficiency by the subse- 
quent fast adiabatic 7t-flipper. Together with the final projec- 
tive measurement along the positive z-direction this gives a 
complete set of measurements of the ±x, ±.y and ±z polar- 
ization components. In this way the final spin state is char- 
acterized with an efficiency close to 100%. The initial de- 
gree of polarization sq = |i(0)| is typically 75%. During the 
evolution so is reduced mainly due to static field inhomo- 
geneities across the storage volume - even without tempo- 
ral fluctuations, i. e. for K{t) = jl^j. Local variations in 
the B-field magnitude cause variations in the Zeeman-energy 
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FIG. 2: Pulses for the B-field (x— , y—, and z— direction) for measur- 
ing the geometric phase: A 7t/2 pulse generates a superposition state. 
This is followed by a cyclic evolution with constant B-field magni- 
tude |6(0I- A subsequent Jt-pulse flips the spin such that the pre- 
ceding unitary evolution can be compensated (By[(])d]). The resulting 
spin state is determined by the polarization analysis (PA). Measuring 
the geometric phase involves a change of the rotation direction 
while keeping the magnitude fixed (By[(|)g]). Identical fluctuations 
are generated in x-direction (K) for measuring geometric dephasing. 



splitting. Consequently, the relative phases in the final spin 
superposition states of the individual neutrons are randomly 
distributed, which leads - on average - to dephasing. This de- 
phasing mechanism causes a further exponential polarization 
loss, s{T) = sqc^^I^^. 72 847(40) ms has been measured 
by a polarization analysis after free precession of the spin su- 
perposition state in a 10 /jT magnetic field, which sufficiently 
exceeds typical evolution times of 500 ms. 

To measure the geometric phase — —£2/2 for K{t) — 
the magnetic guide field pointing initially in the nega- 
tive z-direction is rotated about the x-axis, i. e. By{t) = 
-Bj(0)sin((Of) and B, = B,-(0)cos((Bf) with constant \B{t)\ 
(see Fig.|2|l. An additional offset field Bx in x-direction gen- 
erates a conical section traced out by the magnetic field vec- 
tor and - in the adiabatic limit - also by the spin polariza- 
tion vector. The enclosed solid angle £2 = 7i(l — cosiJ) is de- 
termined by the cone angle = tan^^ B^/B^- To eliminate 
contributions from the dynamical phase (|)^/ we invoke a spin- 
echo scheme 3^. The according evolution path of the 

Tr[d|i+(f))(i+(r)|] of the super- 



spin-up component s^{t) 
position state on the Bloch-sphere is illustrated in Fig. [3 A). 
Depending on the rotation direction the solid angle enclosed 
by the path on the lower hemisphere CIse = ±£i. Thus, if 
the direction of rotation is reversed after a Jt-pulse and the 
field amphtude is kept constant, the geometric phase doubles 
- due to its dependence on the directed solid angle - while the 
dynamical phase cancels. Both the accumulation of the geo- 
metric phase and the cancellation of the dynamical phase has 
been measured using the pulse sequence drawn in Fig. |2]for 
T = 200 ms. The polar angle •& and consequently the solid 
angle £2 is varied by choosing different B^ offset fields. The 
ratio Co/cOi = 2n/{T(0i) w 0.017 ensures adiabaticity of the 
evolution. In Fig. OB) the measured geometric phase <^g is 
plotted as a function of the solid angle £2. The fit to the mea- 
sured data yields (|)^ = —0.5 1 ( 1 )£2 which is in good agreement 
with the expected (|)^ = —£2/2. Residual dynamical phase con- 
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FIG. 3: (A) Path traced out by the spin-up state ?^(/) on the Bloch- 
Sphere while following the adiabatic changes in the magnetic field. 
After the first cycle (path 1-2-3) the state vector is flipped (.v^ I^^) 
and traces out the path 4-5-6 on the lower hemisphere for the second, 
echo cycle. The enclosed solid angle Q. determines the geometric 
phase. (B) Measured geometric phase (filled circles). If the sense 
of rotation is not reversed in the echo pulse, all accumulated phases 
cancel apart from remaining dynamical contributions (solid rect- 
angles). Ideally this yields the same phase as if there was no B-field 
evolution at all ((|)o - open circles). 



tributions, which are not compensated by the spin echo, are 
measured to be 0.22 rad. These are determined by the phase 
difference in the final polarization between the spin-echo with 
identical evolution (Fig.|2]for By[(|)f/]) and without evolution at 
all (B,(f) = const., B,(f) = B,.(f ) = 0). 

For testing the stability of the geometric phase we generate 
field fluctuations K{t) (c. f. Eq. ([T]|) with Lorentzian power 
spectrum 13111 . a bandwidth of F = 100 rads^' and intensity 
a|, = 4 /jT^ as the mean square deviation of the fluctuations. 
We applied a smooth rectangular-shaped window function to 
the noise in the time-domain to avoid non-adiabatic and non- 
cyclic effects. To test the time dependence of the variance 
of the geometric phase <^^^{T) given by Eq. (O the evolution 
time T is changed from T = 35 ms to T = 250 ms. The differ- 
ent fluctuations in subsequent storage cycles lead to different 
phases isf of the final state. But since the noise is identical 
for both first and second part of the spin-echo, these differ- 
ence can only originate in the geometric phase. The dynam- 
ical phase cancels as before in the fluctuation-free measure- 
ments. The average over several storage cycles, i. e. sev- 
eral noise patterns, gives a further shrinking of the length of 
the polarization vector s{T) of the final state additional to the 
unavoidable polarization losses discussed above: In fact, for 



Gaussian distributed (|) we obtain (cos(|)) — exp[— o|/2] cos((|)) 
and same for (sin(|)) in Eq. (|2]). Consequently, the purely geo- 
metric dephasing gives i„(r) = i(r)exp — (40(|)g(r))^ /2 
where the factor 4 is due to the particular type of measure- 
ment f3?l and small fluctuations in are neglected. To sep- 
arate the unavoidable polarization losses from the geometric 
dephasing the geometric phase has been measured with and 
without fluctuations and the ratio of the coiTesponding degrees 
of polarization v,ei(r) = s„{T)/s{T) = exp[- [^(Sifg{T)f /l] 
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gives the variance of the geometric phase a^^. 300 differ- 
ent noise reaHzations have been performed for each value of 
T, where a sequence of six storage cycles forms the polar- 
ization analysis. In Fig. |4] we have plotted the decrease of 
0| ^ as a function of the evolution time T and fixed noise-free 
geometric phase (|)^ = — 2.56rad. The inset shows the corre- 
sponding increase in the relative degree of polarization Vrei- 
Error bars stem from the limited number of noise realizations. 
The solid line indicates the theoretical predictions given by 
Eq. (O for the adjusted experimental parameters without free 
parameters. Due to the low-pass filtering of the coils and in- 
fluences from a non-adiabatic manipulation the data-point at 
r = 35 ms deviates from the theory curve by 3 a. We have 
also verified that the mean geomtric phase remains unaffected: 
A = I (^g) — (j)^ I = ± 0. 1 rad, where denotes the average of 
the geometric phase averaged over the different values of T. 



.07 .035 .023 .017 .014 




50 100 150 200 250 

T [ms] 



FIG. 4: The variance ©^^(r) of the geometric phase as a function 
of the evolution time T in a fluctuating magnetic field. The variance 
decreases for longer evolution times following closely the theoretical 
prediction in Eq. ^ (solid line). The inset shows the increase of the 
degree of polarization Vj-ei relative to the noise-free evolution along 
with the theoretical predictions. 

In summary, we have measured the stability of the adiabatic 
geometric phase with respect to fluctuations in the parameters 
of the driving Hamiltonian as a function of evolution time. A 
spin-echo technique allowed for the observation of the purely 
geometric part of the dephasing of the quantum state by elimi- 
nating dynamical contributions. The acquired data shows very 
good agreement with theoretical predictions and demonstrate 
the vanishing influence of geometric dephasing for slow evo- 
lutions. Clearly, when considering quantum gates a compro- 
mise has to be found between the superior noise resilience but 
slower execution speed compared to dynamical phase gates. 
In this context generalized settings involving non-adiabatic 
geometric phases provide an interesting perspective for 
future experimental efforts. In the adiabatic regime the results 
presented above demonstrate that the geometric phase can in- 
deed be useful for high-fidelity quantum state manipulations. 
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